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ABSTRACT 

We present new observations of the nuclear star cluster in the central parsec of the Galaxy 'with the 
adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our -work allo"ws the 
spectroscopic detection of early and late type stars to niK > 16, more than 2 magnitudes deeper than 
our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We 
find that most of these Wolf Rayet ( WR) , O- and B- stars reside in t'wo strongly 'warped disks bet'ween 
0.8" and 12" from SgrA*, as "well as a central compact concentration (the S-star cluster) centered on 
SgrA*. The later type B stars {rriK > 15) in the radial interval bet'ween 0.8" and 12" seem to be in 
a more isotropic distribution outside the disks. The observed dearth of late type stars in the central 
fe'w arcseconds is puzzling, even 'when allo'wing for stellar collisions. The stellar mass function of the 
disk stars is extremely top heavy 'with a best fit po'wer la'w of d7V/dm oc m~°'^^^°''^. Since at least the 
WR/O-stars 'were formed in situ in a single star formation event ~6 Myrs ago, this mass function 
probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and 
of the early- type stars at distances beyond 12" are compatible 'with a standard Salpeter/Kroupa IMF 
(best fit po'wer la'w of dN/dm oc m~^'^^^°'^). 

Subject headings: Galaxy: center — stars: early-type — stars: luminosity function, mass function 



1. INTRODUCTION 



The central parsec of the Galaxy harbo rs more than 
one hu ndred youn g massive stars (iForrest et al.l 
19871: lAllen et all IToOfil: iKrabbe et al.l [1991|; 
Naiarro et all 119941: iKrabbe et all Il995t iBlum et all 



19951: iTamb lvn et al.l ~IT996[ 



Genzel et all [2 003: C hez et al.l 



Naiarro et al.l 119971: 
20031: iPaumard et all 



20061 : iMartins et al.l |2007| ). This is highly surprising 
since the tidal forces from the central four million solar 
mass black hole should make formation of stars by 
gravitational collapse from a cold i nterstellar cloud very 
difficult if not impossible (jMorrid fl993) . Observations 
have established that most of the WR-stars and O- 
stars (d"warfs, giants, and supergiants) dominating the 
luminosity of the early-type population 'were formed 
in a 'well defined single event ~6 Myrs ago, perhaps 
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as the result of the infall of a g as cloud foUo'wed by 
an in-situ star form ation event (jKrabbe et al.l 119951 : 
IPaumard et"aHl2006f l. About half of these WR/0 -stars 
between 0.8" and 12" from SgrA* appear to reside 
in a 'well defined but highl y 'warped disk that r otates 
clock'wise on the skv (.Levin fc Beloborodo'vl 



Genzel et al 
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20031 IPaumard et al.N2006l: iLu et all 



2003 



2009 



■20091). It is less clear h o'w t he other half 
is dist ributed. IPaumard et all ()2006f ) and iBartko et all 
find evidence for a counter-clock"wise structure. 



perhaps a second disk in a dissolving state. iLu et ahl 
(2009) confirm the first stellar disk but do not observe 
a significant number of stars in the second disk. Within 
about 1" of SgrA* there is a sharp cutoff in the density 
of WR/O-stars. Instead there is a concentration of 
fainter stars 'with randomly oriente d and eccentric 
orbit s : the so-called 'S-st ar cluster' (jEisenhauer et all 
120051: iGillessen et all [20091) . The brighter members of 
this central cusp ar e spectroscopically identified as main 
sequence B stars (i Ghez et al.' 2003; Eisenhaucr et al.1 
[2005; Mar tins et al. 2008; Gillesscn et al. 2009). It is 
currently debated "whether the S-stars formed in the 
disk(s) and subsequently migrated to the central arc- 
second ()Levinl[2007t iLockmann et alll2008f l. or -whether 
they formed outside the central fe'w parsecs, 'were 
injected into near-parabolic orbits by massive perturbers 
and were t hen captured by t he massive black hole (e.g. 
[HiiIi[l988l : rPerets et al.ll200l . 

If star formation indeed has taken place in the vicinity 
of the massive black hole, it is of great interest to ex- 
plore the stellar mass function and spatial distribution 
to gain a better understanding of the processes involved 
in overcoming the tidal barrier. Likewise it is important 
to observationally constrain the characteristics and origin 
of the central S-star cluster. In the following we present 
the results of an extensive imaging spectroscopy survey of 
the central parsec with the adaptive optics (AO) assisted 
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Fig. 1. — The black circles indicate the sample of early-type 
stars with 0.8" < R < 25". The size of the circle scales with 
the K-band magnitude. The large black circles show projected 
distances of 0.8", 3.5", 7" and 12" from Sgr A*. The col- 
ored squares indicate the exposed fields with SINFONI in the 
25 mas/pixel and 100 mas/pixel scale. The color scale indi- 
cates the combined photometric and spectroscopic complete- 
ness for an tuk ~ 14 main sequence early-type star m the 
observed field. The whole inner ~ 12" region is contained in 
lower resolution (250 mas /pixel scale) SINFONI observations 
(Paumard et al\\200di) . The asterisk indicates the visible AO 
guide star. 

integral field spectrograph SINFONI on the ESO/VLT. 
We report new observations aimed at defining the spa- 
tial distribution and dynamics of the nuclear stars up to 
"rriK ^ 16 throughout this region and combine this new 
work with our earlier results on the WR/O-stars and the 
central S-star cluster. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of the imaging survey were carri ed out 
from 2005 to 2009 with SINFONI (EiscnhauereLalJ 
I2003f) on the ESO/VLT, partly using natural guide 
star adaptive optics (AO) an d partly exp l oiting 
the laser guide star fa cility (jRabien et al.l l2003t 
iBonaccini Calia et al.|[2006l ). We covered a total surface 
area of 15 square arcseconds at projected distances below 
4" from Sgr A* with the AO scale (12.5 x 25 mas/pixel 
scale) resulting in a final K-band full width at half max- 
imum (FWHM) of typically about 100 mas. In addition, 
we observed a total surface area of more than 500 square 
arcseconds with the 50 x 100 mas/pixel scale resulting in 
typical K-band FWHMs of about 200 mas. The location 
of the observed fields (typical integration time per field 
2 hours) is indicated in Figure [T] 

We reduced th e data wit h the software SPRED as in 
lEisenhauer et al.l (|l005) and lPaumard et al.l (|2006l ). Our 
primary goal in this project is to determine the spatial 
distribution and dynamics of the early type stars as a 
function of their K-band magnitude / mass. For this 



purpose we need to distinguish from the spectroscopic 
data early type stars from late type stars and to isolate 
the central star cluster from fore- and background stars. 

Late-type stars in the observed magnitude range are 
giants later than KOIII and cool red supergiants. These 
stars show characteristic CO band head absorption fea- 
tures, which can efficiently be used to remove these stars 
from our sample. Only three red supergiant s have been 
obseryed in the central pc of the Galaxy ()Blum et al.l 
I2003t iPaumard et al.l [20060 . most notably IRS 7 at an 
angular distance of 5.5" from SgrA*. In contrast to the 
giants later than KOIII the red supergiants are young and 
massive stars with ages between 7 and a few tens of Myrs 
(|Blum et al.lfl996l) . 

We classified stars as "CO" (late type) or "non CO" 
(early-type) depending on the presence or absence of the 
CO 2 .29/im t; = — 2 band-heads (see also lManess et all 
l2007f ). By adding artificial late type stars to our data we 
determined the average efficiency of the CO detection to 
be better than 98% at tuk = 16.5. In order to rate a "non 
CO" as bona fide early-type star we require the detection 
of at least one of the Br7 , He i 2.058 nm ,or He i 2.112 
/zm lines at a signal to noise level, which allows a radial 
velocity measurement with a maximum uncertainty of 
100 km/s. The sample of bona fide early-type stars is 
essentially free from late type star interlopers. 

Fore- and background early-type stars do not add any 
significant contamination to the observed population of 
early- type stars in the central parsec around SgrA*. 
Due to the large extinction toward the Galactic Cen- 
ter of Ak = 2.8 (see e.g. .Schodel et al.. ,200?! ) . fore- 
ground sources can efficiently be rejected by their very 
blue H-K band colors. There are only two bright in - 
frared detectable foreground stars (jBiretta et al.lll98"^ . 
Background sources are subject to an even higher extinc- 
tion and very red H-K band colors. The cusp of early- 
type stars around SgrA* is very steep (see Figure [2]): 
S(i?) cx For a cluster in equilibrium the distri- 

bution of the line-of-sight distances to SgrA* i s related to 
the s tar surface number density distribution (jAlexanded 
l2005f ). Out of the 136 early- type stars in the interval of 
projected distances 0.8" < R < 12" only 8/3 stars may 
have line-of-sight distances larger than 1/2 pc. 

We corrected the observed K-band magnitude rriK^ohs 
(for the absolu t e magnitude scale calibration see 
iBlum et all 119961 : iManess et all l2007f ) of the stars for 
the variation of extinction over the Galactic Center re- 
gion. We calculated the extinction for each position in 
the Galactic Center field, AK.ohs, by using the observed H 
band and K-band magnitudes of the 20 nearest neighbor- 
ing stars, assuming realistic intrinsic colo rs for the stars 
and a ssuming the extinction law from iRieke fc Riekd 
(fT989l) . A{X) cx A-i-^^ All magnitudes rriK mentioned 
below are referred back to a common standard extinction 
of 2.8: rriK = rriK^ohs + 2.8 - AK,ohs- 

Our total data set contains 177 bona fide early-type 
stars at projected distances up to 25" from SgrA*. Figure 
[Ushows the location of the bona fide early- type stars with 
R > 0.8". Of the 177 bona fide early-type stars in our 
survey 118 are WR/O-stars and 59 are B-dwarfs {niK > 
14.); 28 WR/O-stars and 34 B-dwarfs are reported here 
for the first time). 
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3. NUMBER COUNTS AND COMPLETENESS 
CORRECTIONS 

To study the radial surface density and KLF of the 
early-type stars, we correct the number of observed stars 
for both photometric and spectroscopic incompleteness, 
as well as (for radial surface density plots) for incomplete 
areal coverage. 

The photometric completeness ephot is defined as the 
probability of detecting a star with a given magnitude 
in the image generated by integrating over the spec- 
tral information of the SINFONI data. For this pur- 
pose, we added artificial stars to the spectrally integrated 
SINFONI images, and determined the probability of re- 
detecting the artif icial stars wit h an automatic star de- 
tection algorithm (|Stetsonlll987f ). We used a dense sam- 
pling of the image area with a satisfying number of trials 
(about 2000 trials per image and per flux bin, the exact 
number depending on respective image dimensions) and 
tested the completeness was in flux bins of Am^ — 1. 

The spectroscopic completeness Cspoc is the probability 
of identifying the spectral type of a given, photometri- 
cally detected star in a data cube. For this purpose we 
added and recovered 3D stellar cubes in analogy with the 
photometric completeness estimates. These 3D cubes are 
formed from a spatial PSF with an early-type star spec- 
trum along the wavelength axis. For each photomet- 
rically re-detected star we extracted its s pectrum from 
the d ata cube and computed a CO index (|Maness et al.l 
|2007() . which quantifies the relative depth of the CO band 
head absorption feature. We required a minimum total 
signal to noise ratio of the spectrum in the region of the 
CO band head of at least 5 sigma. If the relative depth of 
the CO band head feature of the artificial early-type star 
spectrum was smaller than twice the noise, the star was 
counted as a "spectroscopically re-identified early-type 
star" . 

We assign errors to both photometric and spectro- 
scopic completeness values. We confirmed that the com- 
pleteness values do not change beyond these errors by 
different settings of the star detection algorithm and the 
CO index calculation. We computed the combined pho- 
tometric and spectroscopic completeness as a function of 
magnitude and location throughout the observed region, 
propagating the individual errors, to the combined com- 
pleteness ecomb(a;, y, Jtik). We calculated completeness 
maps, which give (as a function of magnitude) the com- 
pleteness and its error for each observed point. As an 
example. Figure [T] shows the combined completeness dis- 
tribution for an = 14 main sequence early-type star. 
In our survey an average 50% completeness is achieved 
for an niK — 15 star. The spectroscopy is deeper to- 
ward the central S-star cluster and several fields without 
bright stars (Figure [T]). 

We derived an effective exposed area Aoff (i?i, R^iTtik) 
(including error) over an annulus between i?i and i?2 
around SgrA* as a function of m,K taking the SINFONI 
exposure locations and the completeness values including 
their errors into account: 



Aes{Ri,R2,mK) ^ / ecomb{x,y,mK)dxdy . 

Finally, we computed from the number of observed 
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Fig. 2. — Projected, completeness, and coverage corrected 
surface density as a function of distance to SgrA *. Blue 
points (red triangles) indicate WR/0 (rriK < 14j stars on 
clockwise (counter-clockwise) orbtts and green squares indi- 
cate B-dwarfs in the magnitude bin 14.5 < niK < 15.5. The 
blue and red arrows show the 90% confidence limit to the 
number density of WR/O-stars in the innermost 0.8" and 
counter-clockwise WR/O-stars in the radial interval between 
15" and 25". For comparison, the open asterisks indicate all 
stars with niK < 17 and the full asterisks the identified late 
type stars with tuk < 15.5. 



stars Nsta,rs,ohs.{Ri, R2,''tik) in an annulus between Ri 
and i?2 around SgrA* the K band luminosity function 
KLF{Ri, R2,mK) and the surface density J^{Ri, R2): 



KLF{Ri,R2,mK) 



stars, obs. {Ri,R2,mK) 
Acs{Ri,R2,mK) 



Y^/D D^ ^stars.obs.(-Rl, -R2, 'Tlx) 

i.[Ri, R2) ~ y — - — — — r — 
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4. RESULTS 



4.1. 



Surface density profiles 

Figure [2] shows the projected, completeness, and cov- 
erage corrected surface density profiles for all stars with 
rriK < 17, late type stars with tuk < 15.5, the WR/O- 
stars (rriK < 14), and B-dwarfs in the magnitude in- 
terval 14.5 < rriK < 15.5. The radial surface den- 
sity of the clockwise WR/O-stars follows a power law 
^{R) cx i?-i-4±o.2 between 0.8" and 15". Extrapolating 
this power law to smaller projected distances yields 28±8 
clockwise WR/0 stars for R < 0.8". At the 95% con- 
fidence level we expect more than 12 stars within 0.8". 
We did not observe any WR/0 star in this region. The 
Poisson probability to observe no star for 28 expected 
stars is 5.6 • 10~^^ and the Poisson probability to observe 
none of 12 expected stars is still 6 • 10~^. Therefore we 
conclude that there is a significant cutoff in the distri- 
bution of WR/0 stars at R=0.8". The difference com- 
pared to the R~^ power law suggested by earlier work 
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Fig. 3. — Completeness corrected K-hand luminosity func- 
tions of early-type stars in three radial intervals R < 0.8" 
(red triangles, scaled by a factor 0.05), 0.8" < R < 12" 
(blue points) and 12" < R < 25" (black asterisks). For a 
6 Myr population the best fitting IMF in the radial interval 
0.8" < R < 12", where the disks of early-type stars are most 
prominent, is extremely top heavy and clearly different from 
the IMFs of the S-stars and the field stars beyond 12". It can 
be fitted by a power-l aw IMF with of AN /dm oc m~° ''^*°'^ or 
the IMF proposed bv \Bonnell & RicA 11200^) . The IMF of the 
field stars beyond 12" as well as the S-stars within 0.8" can 
be fitted by a power-law IMF of AN /dm oc m"^'"*^^^"'^, consis- 
tent with a standard Salpeter/Kroupa IMF. These latter KLFs 
can also be fitted by somewhat older populations with a con- 
tinuous star formation history, which would then, however, 
predict lower numbers of mK < 14 stars (relevant for S-star 
cluster). The top x-axis shows the average initial mass of the 
early-type stars contributing to the K band magnitude bins for 
the case of a starburst 6 Myrs ago with an IMF with a slope 



of dN/ dm oc m 



as generated with the STARS code. 



Most of the stars with mK > 14 are B dwarfs, the brighter 
magnitude bins contain O dwarfs, giants and WR stars. The 
main sequence ends at about mK ~ 12, corresponding to an 
initial stellar mass of about 25Mq . 



(|Paumard et al.ll2006t iLu et al.l[2009l : iBartko et alll2009[ ) 

is due to the thorough completeness and coverage cor- 
rection apphed here. The surface density distribution of 
the counter-clockwise WR/O-stars is rather flat between 
0.8" and 15" , mainly because it has a larger central 'hole' 
than the clockwise stars and a maximum in the distri- 
bution at around 4.5" . The radial surface density profile 
of B-dwarfs in the magnitude interval 14.5 < ttik < 15.5 
drops smoothly from the central 0.8", where the S-stars 
reside, out to about 25". Its distribution is similar to 
the clockwise WR/O-stars with a best fitting power law 
of S(i?) oc i?"^-^^"-^. In strong contrast to the early- 
type stars, late type stars with < 15.5 exhibit a flat 
distribution inside of 10", in excelle nt ag reement with 
iBuchholz et"all (|2009[ ) and lDo et all (|2009f ). 

4.2. K Band Luminosity Function 

Figure [3] shows the completeness corrected K-band lu- 
minosity functions of early-type stars in three radial in- 



tervals R < 0.8", 0.8" < R < 12" and 12" < R < 25". 
The slope of the KLF is sensitive to the stellar mass func- 
tion (MF) . If the population formed in situ without being 
affected by dynamical effects that might alter the mix of 
stellar masses in the observed sample, the observed KLF 
of a young population is directly related to the initial 
stellar mass function (IMF). 

The stellar disks (best defined in the radial interval 
0.8" < R < 12", see Figure 2) and the randomly orbiting 
S-stars (in the radial interval R < 0.8") are dynamically 
distinct. The S-stars have a different for mation process 
than the stellar disks (jPerets et al.ll2009[ l. 

The two main scenarios discussed in the liter- 
ature for the origin of the disk-like structurse of 
young sta r s in t he GC are the infallinq cluster model 
j|Gerhardl 120011: iMcMillan fc Portedes ZwartI 
Portcgics Zwa rt etal.1 120031: iKim fc MorriT 



2003 



2003 



Kim et al, i2004l: iGiirkan fc Rasiol 12 005') and the m 
situ star formation model (iLevin fc Bcloborodov 2002 
Genzel et a l. 2003; Goodman 2003; Milosavlicvic fc Loea 
20041 : jNayakshin fc Cuadra 2 005; Paumard et all 
2006; BonncU fc Rice 120081 ; rMapelli et all l2008l ; 
iHobbs fc Navakshin.2009f l. 

In the inf ailing cZwsier scenario, a cluster of young stars 
formed at a distance of at least a few parsecs from the 
GC and then spiraled inward by dyanmical friction. As 
it finally enters the central parsec it is tidally disrupted 
and forms a disk-like configuration. If the cluster is 
dense enough to undergo significant internal mass segre- 
gation during the in-spiral phase this scenario does pre- 
dict a fairly flat MF for the innermost deposited stars 
even if the assumed intrinsic IMF of the cluster is a 



standard one (MF: dN/Ar 



^, IGiirkan fc Rasiol 



120051 ). This is because mass segregation within the 
cluster leads to a differential shedding of stars of dif- 
ferent masses as the cluster is tidally stripped during 
in-spiral. Lower mass stars are deposited further out, 
while the most massive stars sink the farthest in, espe- 
cially if th e cluster contains a cen tral intermediate mass 
black hole (IGiirkan fc Rasioll2005f l. The differential shed- 
ding of different stellar masses combined with the pre- 
dicted flat su rface density of the depos ited cluster stars 
(E cx R-^ '^^. iBerukoff fc Hansenll2006t l would then lead 
to a strong radial change in the MF and a large "sea" 
of B-stars outside the region where most of the O / WR- 
stars are located in the plane of the stellar disks. No 
such gradient is observed in the main part of the disk 
(i? < 12"). Figure 3 does show that the KLF and MF 
of the stars outside 12" steepens, consistent with a nor- 
mal IMF. The observed B-stars at i? > 12", however, 
have large angular offsets to the disk angular momentum 
directions and are compatible with isotropic orbits, indi- 
cating that these stars were not formed in the same way 
as the disk stars. In addition, the steep surface density 
distribution excludes that there are enough B-stars at 
large distances to make up for the shortage of such stars 
in the central 12": 49.3 ± 5.7 stars in the magnitude 
range 13.5 < niK < 16.5 would be needed (assuming 
a standard Salpeter/Kroupa IMF) in the radial interval 
12" < R < 25", but only 13 stars have been observed. 

In the in-situ model, a clump or clumps of gas 
fall into the GC, where they form a disk- like struc- 
ture. The gaseous disk then f ragments to for m stars . 
The detailed observations by iPaumard et al.l (|2006l ): 
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iBartko et all (|2009D : lLii et al.l(|2009[ ) are compatible with 
an in situ formation of these WR/0 stars in a clock- 
wise rotating disk and another highly inclined counter- 
clockwise structure, possibly a disk. The mass segrega- 
tion ti me scales in the GC are much larger than 6 Myr 
([Alexand er .20051) . We therefore assume that the stars 
generated in-situ in the starburst 6 Myr ago are almost 
fully contained in the observed sample of early-type stars 
in the radial interval 0.8" < R< 12". 

Our observations demonstrate clearly that the (I)MF 
in the radial interval 0.8" < R < 12", where the disks of 
early-type stars are most prominent (see Figure [2|), must 
be extremely top heavy. This MF is also clearly different 
from the MFs of the S star cluster and from the region 
outside 12". 

The histograms in Figure [3] show three theoretical 
model luminosity function s with different IMFs, (see e.g. 
iMuench et alll2000l . [200l) . We computed th ese models 
using the population synthe sis code STARS (jSternberd 
[T998t IStCTnberg et all l2008h assuming solar metallicity 
Geneva tracks, for a cluster age of 6 Myr, and an expo- 
nentially decaying star burst with an 1/e time-scale of 1 
Myr. This is the best fitting age and duration of a sin- 
gle star formation event derived from the Hertzsprung- 
Russell diagram distribution of the 0/B-stars a nd the 
ratio s of various sub-types of WR-stars (Paumar d et al.l 
[20061) . STARS computes K magnitudes using empirical 
Schmidt-Kaler bolometric corrections and V-K colors for 
dwarfs, giants, and supergiants. The predicted fractions 
of 0/B main sequence stars and giants as well as WR 
stars agree well with the observations (Paumar d et al.l 
I2006D : Stars with niK < 12 are WR stars and stars with 
rriK > 15 are B main sequence stars. Stars with uik 13 
are WR stars (66%), evolved early Bstars (31%) and late 
Ostars (3%). Stars with ttik ~ 14 a re Ostars (15%) 
and evolved early Bstars (85%). Marti ns et"an ([20071 ) 
presented detailed stellar atmosphere modelings of a sig- 
nificant fraction of the post-main sequence blue super- 
giants and WR stars, which have ages of 4 — 8 Myrs 
and Z AMS masses of at least 6OM0, see also lOtt et al.l 
([19991 ): [Martins et al. (2006). The average properties of 
the generated stars are compatible with the observations. 

The KLF in the radial interval 0.8" < R < 12" can 
be fitted by a power-law IMF of d7V/dm oc rn''°-'^^^°-^ 
ix^ = 4.6 for 6 degrees of freedom). For compari- 
son, the standard Salpeter/Kroupa IMF has a high-mass 
power law of dN/dm oc m~^'^. The IMF found by 
iBonnell fc Ricd ([20081 ) for their IO^Mq cloud simulation 
also gives a reasonable fit to the data (x^ = 7.6, 6 degrees 
of freedom). We have also selected only those early- type 
stars, which have an angular offset of less than 10 de- 
gree s from the lo cal disk angular momentum direction 
(see iBartko et al.l [2009i) and subsequently fitted by the 
theoretical KLF as a function of IMF slope. The ob- 
tained IMF slope agrees with the IMF slope of all stars 
in the radial interval 0.8" < R < 12" within errors. 

The IMF of the field early- type stars beyond 12" and 
the S-stars within 0.8" can be fitted by a power-law IMF 
with a slope of dN/dm oc m~^'^^^°'^, compatible with a 
standard Salpeter/Kroupa IMF. These latter KLFs can 
also be fitted by Salpeter/Kroupa IMFs and continuous 
star formation histories with moderate ages (< 60 Myrs). 

It is straightforward to compute a one-to-one relation 
between observed K band magnitude and initial mass for 




20 40 60 80 100 120 
min angular distance to clockwise disk 

Fig. 4. — Distribution of the reconstructed angular differ- 
ences from the local average angular momentum direction of 
the clockwise system for the early-type stars with projected 
distances between 0.8" and 12". The three distributions cor- 
respond to different K-band luminosities: black: tuk < 14, 
red 14 < rriK < 15 and green mK > 15. Only one out of 11 
stars fainter than rriK = 15 is compatible with the clockwise 
disk. 



stars on the main sequence. However, for evolved stars 
a range of masses contribute at a given K magnitude. 
Therefore we do not estimate the initial mass for each 
star but rather compare the observed distribution of K 
band magnitudes to the simulated magnitude distribu- 
tion for different IMFs. The top x-axis of Figure [3| shows 
the average initial mass of the early-type stars contribut- 
ing to the K band magnitude bins for a 6 Myr old cluster 



for an IMF of dN/dr. 



-0.45 



as generated with the 



STARS code. Most of the stars with tuk > 14 are B 
dwarfs, the brighter magnitude bins contain O dwarfs, 
giants and WR stars. The main sequence ends at about 
rriK — 12, corresponding to stars with initial masses of 
about 25Mo with main-sequence lifetimes of 6 Myr. 

An rriK = 16.5 early-type star for a 6 Myr old pop- 
ulation corresponds approximately to a B5V main se- 
quence star with a ZAMS mass of about 7Mq (see 
Figure [3|) . The most massive stars with an individual 
mass estimates have ZAM S masses of at least 60Afo 
(Martins et al. 200fl, l2007f ). Our estimated IMF slope 
of the stellar disks {dN/dm cx 77i,-n-45±o.3-j therefore at 
least valid over the mass interval 7 — 6QMq. 

4.3. Orbital Angular Momentum Directions 

Figure [H shows the distributions of the reconstructed 
angular differences from the local average angular mo- 
ment um dire c tion o f the clockwise system, as defined in 
Bartk o~et al.l ([2009f ). for the early- type stars with pro- 
jected distances between 0.8" and 12". The three dis- 
tributions correspond to different K-band luminosities: 
mx < 14, 14 < mii < 15, and mx > 15. The peak at 
small angles corresponds to the clockwise system: 50 out 
of 110 stars with mx < 14 have offsets below 20° and 8 
out of 15 stars with 14 < mx < 15 have offsets below 
20°. In contrast, only one out of 11 stars fainter than 
mK = 15 is compatible with the clockwise disk. This 
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Fig. 5. — Cylindrical equal area projection of the distribu- 
tion of significance in the sky (25° aperture, see Bartko et 
al. 2008) for 82 bona fide early-type stars stars (tuk < 14j 
with projected distances between 3.5" and 15". The disk po- 
sitions of Paumard et al. (2006) are marked with full black 
circles. There are two extended excesses visible for clockwise 
and counter-clockwise stars, whith maximum significances of 
8.2 and 7.1 sigma, respectively. We attribute these two ex- 
tended excess systems as two warped systems of stars. 

may be an indication that most of the few later B-dwarfs 
observed in the region of the disks do not belong to the 
clockwise system but rather to the background popula- 
tion seen at projected distances beyond 12", which has 
a more standard IMF. 

According to Figure [2] the clockwise system of mas- 
sive early-type stars is best defined in the radial range 
between 0.8" and 10" and has a steep surface density dis- 
tribution. However, most of the counter-clockwise stars 
reside in the range between 3.5" and 15" , where we have 
added consid erable SINFO NI sky coverage compared to 
the study of iBartko et al.l (2009). In order to test for 
the existence of a significant system of counter-clockwise 
stars, we have computed a sky significance map of angu- 
lar momentum overdensities for all bona fide early-type 
stars stars {rriK < 14) with projected distances between 
3.5 " and 15", se e Figu re [5] For the technical details 
see IBartko et al.l ()2009f ). There are two extended ex- 
cesses visible for clockwise and counter-clockwise stars, 
with maximum significances of 8.2 and 7.1 sigma, re- 
spectively. These maximum excesses are clearly offset 
from the positio ns found at small projected dist ances 
(|Paumard et al.ll2006t iLu et al.|[20"09l : IBartko et"aI]|2009D 
indicating that the two extended excess systems are most 
probably two warped systems of stars. We confirm a sig- 
nificant counter-clockwise system of stars. 

5. DISCUSSION 

The observations discussed in the last section, in 
combination wit h the results pre s ented in recent 
published work fGenzel et al.' 2003; 'Eisenhauer et al] 
2005; Nayakshin & Sunvaev 2005; Schodcl ct al. 2003 
Paumard et al 2006; Lu ct al. 2009; Bartko et al. 20oi 
Gillessen etlH ,2009; .Buchholz et al. 2009 ; Do et all 
20091) vield the following complex and rather unexpected 



properties of the stellar cluster surrounding the central 
black hole, 

• most of the early-type stars in the central parsec 
reside in a combination of a central concentration 
of main-sequence B-stars centered on the black hole 
(the 'S-star cluster'), plus two strongly warped and 
almost orthogonal disks (or planar sets of stream- 
ers) of WR-, O- and early B-stars. One disk is more 
massive than the second (~ lO'' and ~ 5 • 10^ Mg, 
computed from the completeness and exposure cor- 
rected number counts and the best fit IMF). The 
second disk is less well defined, perhaps because 
it is more disrupted. The surface density distri- 
bution of the B-stars appears to decline smoothly 
with approximately constant slope (1.5) from the 
central S-star cluster to the region of the disks. 
However, both the dynamics and the mass func- 
tion change abruptly at 1" . Inside this radius, the 
stellar orbit distribution is random and the ellip- 
ticity distribution is close to but somewhat hotte r 
than a thermal distribution (Gilles sen et al.ll2009r ). 
The observed KLF can be well fitted by a 'stan- 
dard' Salpeter/Kroupa IMF for a single age or con- 
tinuously star forming population of age a few to 
60 Myrs. Just outside of this region most of the 
WR/O-stars are in the dominant clockwise disk, 
which has a 10 degree opening angle and stellar 
orbits with modest (e ~ 0.35) eccentricities. The 
IMF is extremely flat (power law slope 0.45) with 
a mean stellar mass of 30 Mq . Outside of the re- 
gion dominated by the two disks, at R> 12", the 
KLF appears to again approach that of a 'standard' 
IMF. Note, however, that even between 1 and 12" 
~ 20 - 30% of the WR/O-stars, and probably a 
larger fraction of the B-stars do not appear to be- 
long to either disk system. 

• the late type stars of all luminosities, from the 
brightest thermally pulsing AGB stars to the 
faintest, low mass red clump stars, do not show a 
central concentration that might resemble a clas- 
sical equilibrium stellar 'cusp' as predicted by 
all theoretica l studi es (BahcaU & WoH 1976, 197|; 
iFreitag et all 120061; iHopman fc Alexander. .2006i) . 
The brightest late type st ars show a central 
hole of radius a few to 7" (ISellgren et al.l 119901 ; 
IGenzel"eFan[T996l;rHaller et aLlliggl " and even the 
faint early K - giants with ruK ^15 show a flat cor e 
of radius 10" (Suchholz et al.ll2009l; iDo et al.ll2009D . 
see Figure [5] The combination of the S-star cluster 
(mostly unrelaxed early-type stars) and the flat late 
type star distribution 'conspire' to give a shallow 
overall cusp (surface density slope 0. 2... 0.4) found 
in the faint star counts i n Figure [1] (|Genzel et al.l 
[200llSch5del et al.l[200l . 

The observations suggest that the two WR/O/B- 
star disks between 0.8" and 12" are structurally 
and dynamically distinct from both the S-star clus- 
ter and the outer region. The properties of these 
warped disks (or sets of streamers), including in 
particular their steep surface density distribution, 
their top heavy IMF, and their relatively low to- 
tal stellar masses all strongly favor an in situ star 
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forma.tion model j Paum a rd et all l2006l: iLu et al] 
120091 : iBartko et al.li200 9: Bcru koff fc Ha nsen 200"^ 
and broadly agree with the findings in recent hy- 
drodynamical simulations (jBonnell fc Ricd 120081 : 
iHobbs fc NayakshinI l2009h of star formation trig- 
gered by gas cloud infall. The sharp transition at 
0.8" between disk zone and S-star cluster, both in 
terms of dynamics and mass function, in our opin- 
ion also strongly disfavors migration scenarios from 
the disks for the origin of the S-star cluster. In 
turn this sharp transition supports the 'injection 
and capture' scenarios mentioned in the Introduc- 
tion. 

The large central core (and perhaps even cen tral 
hole, see iBuchholz et "all 120091 : iDo et al.ll2009t ) in 
the density distribution of late type stars of all 
luminosities is puzzling and currently not under- 
stood. There are a number of possible interpreta- 
tions, none of which at prese nt offers a compe lling 
explanation. A s discu ssed by iDo et al.1 (|2009D and 
IBuchholz et al.l ()2009t ) equilibrium mass segrega- 
tion by itself canno t account for this distr i butio n 
(jFreitag et al.l 120061: iHopman fc Alexan der '2006D. 
as even in a multi-mass cluster the lightest stars 
attain a radial density distribution that is steeper 
than (probably excluded by the data: iDo et al.l 
|2009[ ). Simulations indicate that physical colli- 
sions with main-sequence stars and stellar black 
holes can remove moderately bright giants in the 
central arcsecond, but n ot over a rnuch la rger re- 
gion, nor to rriK ~ 15 (jPale et al.l l2009f ). Like- 
wise, tidal disruption of s tars may play a role 
near the massive black hole ()Perets et al.ll2009t ) but 
it needs to be investigated whether they are fre- 
quent enough to remove the entire old cusp. The 
in-spiral of an intermediate mass black hole can 
plausibly gouge ou t a large enough core in the 
stellar distribution (jMilosavlievic fc MerrifrB 120011 : 



iBaumgardt et al.l l2006f) but there is no evidence 
for such a second black hole from the dynamical 
properties of the S-stars or late type star s, nor 
from t he motion of SgrA* it self fGillcsscn et all 
20091: [Gualandris & Mcrritt 2009; Triooe et aj 
20081 : iReid fc Brunnthaler 2004). M errit t (20091) 



has pointed out that the two-body relaxation time 
scale probably is significantly longer than the age 
of the Galactic Center star cluster (and the Hub- 
ble time) throughout the central parsec, especially 
if the spatial distribution of the giants is similar to 
that of most of the stellar mass. As a result the ob- 
served large core may reflect the initi al conditions 
of the Galactic Center nuclear cluster. iTrippe et al.l 
(|2008D find, however, that the dynamics of the 
old star cluster is consistent wit h a relaxed, fully 
phase mixed system. In addition I Alexand"ei1 ()2007f ) 
has pointed out the importance of relaxation pro- 
cesses much faster than the standard two-body 
rate (inclu ding the massive perturbers proposed by 
iPerets et a l. 2007). The top-heavy IMF discussed 
in this paper (if applicable to early star formation 
episodes) would obviously also lead to a lack of 
old low mass giants in the core. To expla i n the 
flat or even inverted radial slope (|Do et al.l 120091 : 
[l3uchholz ct al. 2009,, this paper), the IMF proba- 
bly would have to depend strongly on radius. The 
remarkable properties of the Galactic Center nu- 
clear star cluster remain puzzling and continue to 
give us food for thought and further work. 
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